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Abstract. The first complete far- infrared spectrum of T 
Tau has been obtained with the LWS spectrometer on- 
board the Infrared Space Observatory (ISO), which de- 
tected strong emission from high-J (J=14-25) CO, para- 
and ortho-H20 and OH transitions over the wavelength 
range from 40 to 190 /im. In addition the [OI]63^m, 
[OI]145/im and [CII]158/im atomic lines were also de- 
tected. Most of the observed molecular emission can be 
explained by a single emission region at T'--^300-900 K 
and ~10^~^cm~'^, with a diameter of about 2-3 arc- 
sec. This corresponds to a very compact region of 300 - 
400 AU at the distance of 140 pc. An higher temperature 
component seems to be needed to explain the highest exci- 
tation CO and H2O lines. We derive a water abundance of 
1-7 -10^^ and an OH abundance of ^ 3-10^^ with respect 
to molecular hydrogen, implying H2O and OH enhance- 
ments by more than a factor of 10 with respect to the 
expected ambient gas abundance. 

The observed cooling in the various species amounts 
to 0.04 L0, comparable to the mechanical luminosity of 
the outflow, indicating that the stellar winds could be re- 
sponsible of the line excitation through shocks. 

In order to explain the observed molecular cooling in 
T Tau in terms of C-type shock models, we hypothesise 
that the strong far-ultraviolet radiation field photodisso- 
ciates water in favour of OH. This would explain the large 
overabundance of OH observed. 
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The estimated relatively high density and compact- 
ness of the observed emission suggest that it originates 
from the shocks taking place at the base of the molecu- 
lar outfiow emission, in the region where the action of the 
stellar winds from the two stars of the binary system is 
important. 

Key words: Stars: formation; Stars: individual: T Tau; 
Stars: pre-main sequence; ISM: individual objects: T Tau; 
ISM: jets and outfiows; Infrared: ISM: lines and bands 



1. Introduction 

T Tau has been extensively studied from radio to ultravi- 
olet wavelengths because it has long been considered the 
prototype of a class of pre-main sequence stars. In recent 
years it has become clear that T Tau is in reality a very 
complex system and that it differs from other T Tauri 
stars. It is in fact known to be a binary system (Dyck et 
al. 1982), containing an optical stellar component (T Tau 
N) and an infrared companion 0.7 arcsec to the south (T 
Tau S), corresponding to 100 AU at the 140 pc distance of 
the Taurus Auriga dark cloud (Wichmann et al. 1998). T 
Tau S dominates the bolometric luminosity of the system 
(Chez et al. 1991), but it has no detectable optical coun- 
terpart to a limiting magnitude of V=19.6, suggesting an 
optical extinction greater than 7 mag (Stapelfeldt et al. 
1998). 

Since Herbig's (1950) optical spectra of Burnham's 
nebula, surrounding T Tau several arcsecond across, re- 
sembling a so-called HH object, it was clear that the in- 
teraction of stellar winds with the surrounding molecular 
medium is at work. Molecular outflow activity was first 
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mapped by Edwards & SncU (1982) in ^^CO J=l-0 and 
J=2-l, who found that 95% of the high velocity molec- 
ular gas is associated with blueshifted material. The di- 
rection of the detected outflow is roughly parallel to the 
line of sight, but the emission also shows a region extend- 
ing 2 arcmin to the south and east of T Tau N with a 
secondary peak in the blueshifted wing. Higher resolution 
maps of the ^^CO J=3-2, J=6-5 and C^^O J=l-0, J=2-l 
and HCO-I- emission later showed (Schuster et al. 1993, 
Momosc et al. 1996, Schuster et al. 1997, Hogerheijde et 
al. 1998) a complex outflow system that could originate 
from the different components of the binary system. 

The fast stellar winds observed through forbidden opti- 
cal line emission (Bohm & Solf 1994), revealed five distinct 
kinematic components that suggest that both the primary 
star and the companion may drive separate bipolar out- 
flows. A giant Herbig-Haro flow was recently discovered 
(Reipurth et al. 1997) around T Tau and is interpreted as 
originating several thousand years ago from T Tau S. 

Strong and extended H2 ro-vibrational emission was 
found quite early around T Tau (Beckwith et al. 1978). 
Recently, high resolution H2 imaging (Herbst et al. 1996, 
1997) indicated that the extended molecular hydrogen 
emission arises from the impact on the ambient cloud of 
two outflow systems oriented NW-SE and E-W. These 
originate from the two stars, each with its circumstellar 
disk, and the emission is distributed equally over T Tau N 
and T Tau S. Infrared adaptive optics observations in H2 
show instead that the emission is concentrated on T Tau S 
and is interpreted in terms of shocks occurring as matter 
accretes onto the circumstellar disk of T Tau S (Quirren- 
bach & Zinnecker 1997). No firm conclusion is therefore 
reached on this problem. 

T Tau has associated a substantial amount of mass of 
dust, it is therefore luminous in the millimeter continuum 
(Adams et al. 1990, Beckwith et al. 1990). A circumstellar 
disk has been detected both with CO interferometry and 
infrared scattered light (Weintraub et al. 1989, Momose 
et al. 1996). Later, millimeter continuum interferometry 
at 0.9 and 3mm (Hogerheijde et al. 1997, Akeson et al. 
1998, respectively) was used to derive a total mass of 0.04 
Mq for the circumstellar disk around T Tau N and at 
least 10 times smaller for that associated to T Tau S. A 
circumbinary envelope would also be required to fit the 
continuum energy distribution. 

Far-infrared spectroscopy provides powerful diagnostic 
lines from abundant molecular species like CO, H2O and 
OH, that can be used to clarify the physical processes 
at work in the complex T Tau system. In this paper we 
present the far-infrared spectrum observed from the Long 
Wavelength Spectrometer (LWS, Clegg et al. 1996) on- 
board the Infrared Space Observatory (ISO, Kessler et al. 
1996). Additional data from the Short Wavelength Spec- 
trometer (SWS, de Graauw et al. 1996) are also used for 
discussing the molecular emission properties. The main 
results of the SWS are presented by van den Ancker et 
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Fig. 1. The complete ISO LWS and SWS spectrum of T 
Tau. The IRAS data are also shown for comparison 

al.(1999). The continuum far-infrared spectrum of T Tau 
will be discussed in a forthcoming paper (Pezzuto et al. in 
preparation). 

2. Observations 

T Tau has been observed with ISO (Infrared Space Ob- 
servatory) using the LWS (Long Wavelength Spectrome- 
ter, Clegg et al. 1996). A full low resolution (R ~ 200) 
spectrum of the source from 45 to 197 fjbm was obtained 
during revolution 680, corresponding to September, 25, 
1997. The beamsize is on average 80 arcsec, depending 
on the wavelength. The spectrum was made up of 23 full 
grating scans oversampled at 1/4 of a resolution element, 
with each spectral sample integrated for 11.5 sec, with a 
total integration time of 4265 sec. Besides the observa- 
tion on-source, full grating spectra were also collected at 
four off-source positions. In Table 1 we present the jour- 
nal of the LWS observations, which includes source and 
off-source positions and total observing time (OTT). 

The raw data were reduced and calibrated using ver- 
sion 7 of the LWS pipeline, which achieves an absolute ac- 
curacy of about 30 % (Swinyard et al. 1998). Post-pipeline 
processing was carried out with the ISAP package and in- 
cluded removal of spurious signals due to cosmic ray im- 
pacts and averaging the grating scans of each detector. 

Besides the data measured by the LWS, we also discuss 
in this paper the detection of H2O and OH emission lines 
observed by the SWS. The details of these observations 
are reported by van den Ancker et al.(1999). 
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Fig. 2. The ISO-LWS emission line spectrum of T Tau, from which the continuum has been subtracted. 



3. Results 

The 2-200 /im far-infrared spectrum of T Tau composed 
by the ISO LWS and SWS spectra is shown in Fig. 1. 
The displacement (of ^ 15%) between the flux level of the 
two instruments at 45/im is well within the calibration 
uncertainties. In the figure, IRAS photometry is also re- 
ported for comparison. As discussed in van den Ancker et 
al. (1999), the higher ISO fluxes compared to the IRAS 
data (30-50%) can be explained by the flare that occurred 



to the system in 1990-1991 (Ghez et al. 1991; Kobayashi 
et al. 1994), after which the infrared luminosity of T Tau 
did not return to the pre-outburst value. 

The line spectrum (Fig. 2) is very rich in molecu- 
lar emission lines from the rotational spectra of carbon 
monoxide, water and hydroxyl. The line fluxes, computed 
by fitting gaussian profiles to the lines, are listed in Tables 
2-5. All of the CO transitions with Jup = 14-25 appear in 
the spectrum. However, we cannot assign a flux, but only 
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Table 1. Journal of the ISO-LWS observations of T Tau 



Pos. 


2000.0 Coordinates 


OTT 




R.A. DEC. 






h m s dcg ' " 


(sec) 


T Tau on 


4:21:59.4 +19:32:06.5 


4265 


T Tau off N 


4:21:59.4 +19:33:46.5 


1345 


T Tau off S 


4:21:59.3 +19:30:26.5 


1345 


T Tau off W 


4:22:06.4 +19:32:06.0 


1345 


T Tau off E 


4:21:52.4 +19:32:07.0 


1345 



an upper limit, to the Jup = 23 Une because it is blended 
with the water line 0-H2O 4i4-3o3. Water lines are ob- 
served both in the ortho (o-) and para (p-) form. Most of 
the back-bone lines up to excitation temperatures of more 
than 700 K arc detected for ortho-H20 and up to nearly 
400 K for the para-H20, as well as most of the strongest 
transitions falling in the LWS spectral range. Also the 
strong OH lines arc detected and arc particularly strong, 
with excitation temperatures up to 600 K. [OI]63/xm and 
145/um and [CII]158/im are the only atomic lines present 
in the ISO-LWS spectrum. [OI]63/im emission from T Tau 
was already detected with the Kuiper Airborne Observa- 
tory from Cohen et al. (1988), who measured a total flux 
in a 47" aperture of (242 ± 25) x lO'^o W cm'^, which 
is consistent with our value, which is however much more 
precise. This indicates that the flare occurred in 1990-91 
did not influence the far-infrared line emission. In the four 
off-source positions only the [CII]158/im line was detected 
(see Sect. 3.4 and Table 5). 

A large velocity gradient (LVG) code, solving the level 
population equations in a plane parallel geometry (Nisini 
et al. 1999a, Giannini et al. 1999) was used to model 
the observed molecular emission line intensities from CO, 
H2O and OH. As a first approximation, the local radiation 
field was not taken into account in the radiation transfer 
calculations. We will however discuss in Sect. 3.3 what is 
the effect of considering the infrared local radiation field 
in the OH model. 

The model has many free-parameters (gas tempera- 
ture and density, intrinsic width of the line, column den- 
sity and emitting area or filling factor, which is related to 
both number and column densities) that cannot be eas- 
ily constrained simultaneously. Given the assumption that 
all of the molecular lines originate from the same emitting 
gas, we started our analysis by deriving a range of tem- 
perature and density that is allowed from fitting the CO 
lines, that are likely optically thin. These parameters are 
then used to model the water and OH lines, which, unlike 
the CO lines, have high optical depths due to their strong 
radiative transitions in the far-infrared. 



Table 2. Measured CO line fluxes from the LWS grating 
spectrum, with la uncertainties. Upper limits are at 3(7. 
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(/im) 


TIT 

±u vv 


Clll I 




CO 


28-27 


yo.oo 








CO 


27-26 


yo. 1 / 








CO 


26-25 


100.46 


<5.3 




104.40 


CO 


25-24 


104.44 


7.0 


1.3 


108.81 


CO 


24-23 


108.76 


6.7 


1.1 




CO 


23-22 


113.46 


<10.5* 




118.58t 


CO 


22-21 


118.58 


5.1 


0.8 


124.31 


CO 


21-20 


124.19 


5.2 


0.9 


130.40 


CO 


20-19 


130.37 


7.7 


0.5 


137.17 


CO 


19-18 


137.20 


9.6 


0.7 


144.78t 


CO 


18-17 


144.78 


12.6 


0.5 


153.24 


CO 


17-16 


153.27 


14.1 


0.4 


162.81t 


CO 


16-15 


162.81 


14.0 


0.7 


173.63t 


CO 


15-14 


173.63 


15.5 


1.0 


185.93 


CO 


14-13 


186.00 


14.3 


1.0 



Notes: f: wavelength was fixed for deblending. 

★: this line is blended with the 0-H2O 4i4-3o3 (see text), the 

total flux has a Icr uncertainty of 0.5 • 10-^° W cm-^. 

3.1. CO emission 

For the CO model, we computed the coUisional downward 
rates for levels with Jup < 60 and T > 100 K using the 770 
coefficients taken from McKee et al. (1982), while the up- 
ward rates were computed using the principle of detailed 
balance. Radiative decay rates were taken from Chacke- 
rian & Tipping (1983). 

The distribution of the observed CO line fluxes as a 
function of the rotational quantum number is shown in 
Fig. 3. Because the CO lines are optically thin, their emis- 
sion, in the LVG model considered, does not depend on 
the velocity gradient and thus on the assumed line-width. 
We have considered for our fit only the transitions with 
Jup less than 22. Our data are consistent with gas tem- 
peratures ranging from T = 300 to 900 K and molecular 
hydrogen densities of n = lO^-^cm"^. The two extreme 
models consistent with the data have: 

T = 300 K and uh^ = 4 • lO^cm-^; 
T = 900 K and uh^ = 2 ■ lO^cm-^. 

Fig. 3 shows that the transitions with Jup = 24 and 25 
have a flux level which is too high to be explained by the 
same gas component of the other lines and may indicate 
the presence of a warmer gas emission. This warmer com- 
ponent, which cannot easily be constrained by the higher 
Jup transitions observed, could also affect the Jup = 21 
and 22 lines. However fitting the component covering the 
lines 14 < Jup < 20, results in the same parameters as the 
low temperature model above. 
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Fig. 3. Model fits tlirough tiie observed CO lines. The 
range of temperatures and densities compatible with the 
observations are indicated. The higher observed J lines 
(J=24, 25 and 26) have fluxes too high to be fitted by 
the same parameters as the other lines, suggesting the 
presence of a second component. 

3.2. H2O emission 

As outlined in the previous Sect., we adopted the temper- 
ature and density as derived from the CO lines models to 
fit the observed H2O line fluxes. 

We considered in the computation 45 levels for both 
the ortho and para species (i.e. excitation energies up to 
~ 2000 K): radiative rates are taken from Chandra et al. 
(1984) while the H2O-H2 collision rates are derived from 
Green et al. (1993). We assumed an ortho/para abundance 
ratio of 3, equal to the ratio of the statistical weigths of 
their nuclear spins. 

The other parameters that enter in the model are the 
velocity gradient in the region (dV/dr) and the projected 
area of the emission region. The optical depth in the 
lines is directly proportional to the water column den- 
sity (N(H20)). Since the ratios of different lines depend 
on their relative optical depths, we can use them to con- 
strain dV/N(H20). On the other hand, the absolute line 
intensity depends on both the column density and the 
projected area of the emission region, and therefore if we 
assume a velocity linewidth dV, we can estimate both the 
column density and the emission region size. 

The results of the model fitting are shown in Fig. 4 
and Fig. 5 for ortho- and para-H20 respectively. Almost 
all the lines in the LWS wavelength interval are well re- 
produced by the model. The differences in the flux level 



Fig. 4. Comparison of the modeled ortho-H20 line fluxes 
(filled triangles) with those observed (open circles) for the 
two models considered. 
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Fig. 5. Comparison of the modeled para-H20 line fluxes 
(filled triangles) with those observed (open circles) for the 
two models considered. 

predicted in the two extreme models are quite small, indi- 
cating that the H2O emission is not very sensitive to the 
exact value of temperature and density in the range. On 
the other hand, we note that the ortho-H20 lines short- 
ward of 50/um are not fitted by our models. These lines. 
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Table 3. Measured water line fluxes from the SWS (up- 
per list) and LWS (lower list) grating spectrum, with la 
uncertainties. Upper limits are at 3a. 
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Notes: f: wavelength was fixed for deblending. 

this line is blended with the CO 23-22 line (see text), the 
total flux has a la uncertainty of 0.5 • 10-^° W cm-^. 
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Fig. 6. Comparison of the modeled OH line fluxes (filled 
triangles) with those observed (open circles). 



Using this emission area, the CO column density that 
we derive from the observed CO absolute line fluxes is 
N(CO) = (0.7-2.0)1018 cm-2 and therefore an H2O/CO 
abundance ratio of ~ 0.1 - 0.7. Assuming a standard CO 
abundance of 10"^, the water abundance with respect to 
H2 is ^ (1 - 7) • 10"^. This value implies an enhancement 
with respect to the expected abundance in the ambient 
gas of at least a factor of 10 (e.g. Bergin et al. 1998). 



originated by levels at energies higher than 500 K, are 
brighter than our predictions, indicating that a warmer 
component might be required, as also suggested from the 
higher J CO lines. 

An estimate of the intrinsic lincwidth dV can be given 
if we relate the observed emission with the outflow/ wind 
activity taking place in the close environment of the T 
Tau binary system. The molecular outflow has been traced 
by different lines at near infrared and millimeter wave- 
lengths (H2, CO, HCO+), showing linewidths of a few km 
s~^. Adopting a velocity of 10 km s"-'^, close to the out- 
flow velocities of 7.9 km s-i(red lobe) and 9.7 km s-i(blue 
lobe) measured by Levrault (1988) and those measured by 
Hogcrhcijde et al. (1998) in ^^CO 3-2 (of 12. and 15.6 km 
s-ifor the red and blue lobes , respectively), we derive 
a water column density of (2-5) 10^^ cm"^, while the pro- 
jected area is (4-9) arcsec^. This corresponds to a diameter 
of only 300 - 400 AU, assuming spherical simmetry. The 
compactness of this emission region will enable us to put 
constraints on the physical mechanisms responsible of the 
observed emission (see Sect. 4). 



High H2O abundances are common in young stellar 
objects: the ISO spectrometers have in fact found strong 
emission from gas-phase H2O from massive young stars 
(Harwit et al. 1998, Gonzalez-Alfonso et al. 1998) and 
from low mass outflow driving sources (Liseau et al. 1996, 
Saraceno et al. in preparation, Ceccarelli et al. 1998) with 
abundances in the range 1-5T0~^, rising to values as high 
as ~ 5-10-^ in L1448mm (Nisini et al. 1999b) and Orion 
(Harwit et al. 1998). 

3.3. OH emission 

For the OH models, we have considered 20 levels. The 
collisional downward rates are from Offer & van Dishoeck 
(1992) and the radiative decay rates are from the HITRAN 
catalogue (Rothman et al. 1987). 

If we adopt the same parameters as derived from the 
above analysis also for the OH, we find that a better agree- 
ment between data and models is achieved with the lower 
temperature model (T = 300 K). The estimated OH col- 
umn density is N(OH) ~ 4 • 10 cm^^ and therefore a 
X(OH) ~ 2.7 • 10"^ The results of the OH model fitting 
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Table 4. Measured OH line fluxes from the SWS (upper list) and LWS (lower list) grating spectrum, with la 
uncertainties. Upper limits are at 3c7. 
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1.5 



OH 2ni/23/2-2ni/23/2 53.26/53.35 

OH 2ni/29/2-2ni/27/2 55.89/55.95 17.9 3.3 

OH 2n3/29/2-2n3/27/2 65.13/65.28 23.7 1.0 

OH 2ni/27/2-2ni/25/2 71.17/71.22 18.2 3.3 

OH 2ni/2l/2-2ni/23/2 79.12/79.18 17.8 3.0 

OH 2n3/27/2-2n3/25/2 84.42/84.60 25.9 2.4 

OH 2ni/25/2-2ni/23/2 98.73 9.8 1.0 

OH 2n3/25/2-2n3/23/2 119.23/119.44 14.9 1.5 

OH ^ni/23/2-^ni/2-l/2 163.12/163.40 8.9 0.7 



53.17 


0.75 


55.91t 


0.35 


65.23 


0.29 


71.22 


0.30 


79.16t 


0.29 


84.59 


0.59 


98.73t 


0.6 


119.44 




163.12t 


0.6 



Notes: f: wavelength was fixed for deblending. 



are shown in Fig. 6. As can be seen from this figure, not all 
the lines can be reproduced by the models, in particular 
both the 163/im line and those shortward of 60/zm. 

Because T Tau is relatively bright in the continuum at 
the far infrared wavelengths (sec Fig.l), such a discrep- 
ancy could be due to the pumping from the local ther- 
mal radiation field. To test this possibility we also com- 
puted models for the OH transitions including the field 
originated from dust at a temperature of 300K, using the 
model of Cesaroni & Wamsley (1991). We found that the 
inclusion of the local infrared radiation field increases the 
emission in the lines with A less than 100/xm and the 
163/im line. 

3.4- Atomic emission 

The detection of [CII]158/zm in the four off-source posi- 
tions around T Tau (see Table 5) clearly shows that most 
of the ionized carbon emission (~ 80%) is from an ex- 
tended region and not originated in the vicinity of T Tau. 
The intrinsic emission in the LWS beam centered on T 
Tau is about 2 • 10"^° W cm-^. This implies a ratio 
[OI]63/im/[CII]158/xm of 115, greatly in excess of that 
expected from photodissociation region models (Kaufman 
et al.l999; Burton et al. 1990). 

Finally the ratio of [OI]63/im/145/im=28.3 is such that 
no oxygen self- absorption should occur, as it often appears 
to be the case towards pre-main sequence sources (Sara- 
ceno et al. 1998). This suggests that there is no cold gas 
in front of the source, in accordance with the geometry of 
the outflow directed towards the observer. 

As outlined in van den Ancker et al. (1999), based on a 
larger set of fine-structure lines detected at shorter wave- 
lengths, we argue that the atomic line emission observed is 



consistent with the presence of J-type dissociative shocks. 
On the other hand, the excess [Cn]158/L<m emission on- 
source could also be due to a local photodissociation region 
(PDR) , possibly originated from the far-UV field of T Tau 
(see Sect. 4). 



3.5. Total cooling 

Table 6 summarises the physical quantities derived from 
the observed molecular spectra of CO, H2O and OH, 
adopting the two models considered: the column densi- 
ties and the total cooling luminosities are given for each 
molecular species. We also give the observed values of the 
cooling derived from the sum of all the detected fluxes. For 
deriving Lh2, we used the line fluxes reported in van den 
Ancker et al. (1999) and the line H2 1-0 S(l) at 2.12/im, 
given in Carr (1990). A comparison between the observed 
and modeled cooling shows that the observations of CO, 
H2O and OH can account for most of the modeled cooling. 
The underestimate of the OH cooling by the 900 K model 
confirms that this latter is probably inadeguate to explain 
the observed OH emission. The total radiated cooling ob- 
served from these species, including [01], sums up to about 
0.04 Lq, and has to be considered as a lower limit. Tak- 
ing the outflow parameters from the literature (Levreault 
1988; Mundt 1984) we derive a total mechanical luminos- 
ity of about 0.05 Lq (we have taken an average veloc- 
ity of 10 kms"^, a total outflow mass of 0.22 Mq, and a 
dynamical timescale of 40,000 years). The radiative lumi- 
nosity observed in the far-infrared is therefore comparable 
to the outflow mechanical luminosity. This is expected if 
the stellar winds from the stars are driving the outflows 
and the shocks, traced by the far-infrared lines, accelerate 
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Table 5. Measured atomic line fliixes from the LWS grating spectrum with uncertainties. 



Pos. 


Ao6s 


Line id. 






F 


AF 




(yum) 






(/urn) 


(io-2» 


W cm-2) 


on 


63.21 


[OI] ''Pi ^ 


'P2 


63.18 


230.6 


1.0 


on 


145.52t 


[OI] ^ 

[CII] ^P3/2 ^ 
[CII] 2P3/2 ^ 
[CII] 2P3/2 ^ 


'Pi 


145.52 


8.15 


0.78 


on 


157.76 


'Pl/2 


157.74 


10.6 


0.5 


off-N 


157.80 


'Pl/2 


157.74 


7.9 


0.7 


off-S 


157.79 


^PV2 


157.74 


8.8 


1.2 


off-W 


157.76 


[CII] 'P3/2 ^ 
[CII] 2P3/2 ^ 




157.74 


8.8 


0.9 


off-E 


157.74 


'Pl/2 


157.74 


9.1 


0.7 



Notes: f: wavelength was fixed for deblending. 

the ambient medium into the molecular outflow (Davis & 
EislofFel 1995). 



Table 6. Physical parameters of the molecular and atomic 
emission 





observed* 


"lower T" 


"higher T" 


Temperature T(I< 




300 


900 


Density nnj (cm~ 




4-10^ 


2-10^ 


Nco (cm"^) 




2 • 10^* 


7 • 10" 


NH2O (cm^^) 




2 • 10" 


5 • 10" 


NoH (cm-2) 




4-10" 


4 ■ 10" 


Leo (Lq) 


0.0067 


0.0096 


0.010 


LhsO (L©) 


0.008 


0.014 


0.009 


LoH (L0) 


0.012 


0.019 


0.009 


L[oi] (L©) 


0.014 






Lh2 (Lq) 


0.005 







Notes: The observed cooling is computed by summing all 
the detected lines. 



4. Discussion 

Once the physical conditions of the emitting gas in the 
vicinity of T Tau have been established, we can now pro- 
ceed to compare the results with existing models of shock 
excitation. The far-infrared line emission of young stel- 
lar objects is mainly originated from two physical pro- 
cesses: the excitation from photoionizcd and photodisso- 
ciated (PDR) regions (Tielens & HoUenbach 1985) and 
the shock excitation produced by the interaction of su- 
personic winds with the ambient medium. Depending on 
wind velocity, magnetic field and ion density, two kinds 
of shocks with different far-infrared spectra are predicted 
from models: 

i) high velocity dissociative J shocks (e.g. HoUenbach & 
McKee 1989), in which temperature, density and veloc- 



ity have a discontinuous jump (J) on the shock front, 
molecules are dissociated and atomic lines are the dom- 
inant coolants; 
ii) low velocity non-dissociative C shocks (e.g. Kaufman 
& Neufeld 1996, Draine et al. 1983) in which the ion 
Alfven velocity is larger than the shock velocity and 
the magnetic field transmits energy faster than the 
shock velocity; in this case temperature, density and 
velocity have a continuous (C) variation and molecules 
are the dominant coolants. 

Instead of using intensities of many tens of different 
lines, we can obtain a better comparison between our data 
and the shock model predictions using the total cooling 
from a single species. 

In Fig. 7 we show the water cooling as a function of 
the [OI]63/zm cooling, both normaHzed to the high- J CO 
cooling. The C-type shock models of Kaufman & Neufeld 
(1996) are considered. The J-type shock models (HoUen- 
bach & McKee 1989) are also indicated (see figure caption 
for details). Together with the position of T Tau, we also 
show in this figure the positions of L1448 (Nisini et al. 
1999b), IC1396 (Saraceno et al. 1999), IRAS16293-2422 
(Ceccarelli et al. 1998), and the Herbig-Haro objects HH54 
(Liseau et al. 1996), HH25 and HH26 (Benedettini et al. 
1998). 

T Tau is in a central position, showing that both C- 
type and J-type shock models could explain the obser- 
vations. It has to be noted, however, that according to 
Kaufman & Neufeld models, its position implies a shock 
velocity between 10 and 15 km s~^, in a range where wa- 
ter production is triggered but it is not at its maximum 
efficiency. We can see in the figure that other pre-main 
sequence sources also cluster in the same region of T Tau, 
indicating that these shock conditions are fairly common 
in the environment of young stellar objects (Nisini et al. 
1998). 

Fig. 8 shows the water cooling as a function of the OH 
cooling , both normalized to the high-J CO cooling. As 

before, we consider both J-type and C-type shocks. The 
position of T Tau in this plot appears to be consistent 
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Fig. 7. Total cooling from water lines versus 
[OI]63/Ltm cooling, both normalized to the CO (high- J) 
cooling in few objects, including T Tau (filled triangles) 
and according to shock models. The open triangle shows 
the position that T Tau would have if the water cooling 
were increased by a factor 3. C-type shock models of 
Kaufman & Neufeld (1996) are shown in a grid where 
shock velocity increases from the right to the left from 
10 < Vs < 40 km s~^ (dashed lines) and density from the 
bottom to the top from 10^ to 10^'^ cm^"^ (solid lines). J- 
type shock models (Hollenbach & McKee 1989) are shown 
for comparison as large dashed circles, where S is the 
standard model, with preshock density of Uq = 10^ cm~^, 
shock velocity of Vg = 80 km s~^ and magnetic field of 
B = 158/iG; UV-: the far-ultraviolet field is reduced by a 
factor of 10; B-: the magnetic field is reduced by a factor 
10; B+: the magnetic field is increased by a factor 10; 
GR-: the grain size distribution is extended down to 10 
A; HL: H2 formation on grains is equal to zero; HZ: H2 
formation heating is set to zero. 



with J-type shocks. The main reason why C-shocks fail 
to reproduce the observed cooling is the overabundance of 
OH molecules, by at least a factor 10. Only a model which 
includes the efifects of the presence of a high UV field from 
T Tau itself can reproduce the observed values. 

Fig. 9 shows the water cooling as a function of the H2 
cooling, both normalized to the high-J CO cooling. T Tau 
lies exactly in the center of the region of C-type shocks 
and cannot be explained by J-type shocks. As expected, 
the strong H2 emission cannot be accounted for by J-type 
shocks. Wc therefore rule out the hypothesis that the ma- 
jor responsible of the observed excitation are J-shocks. 



Fig. 8. Total cooling from water lines versus total cooling 
from OH lines, both normalized to the CO (high-J) cool- 
ing in few objects, including T Tau (filled triangles) and 
according to shock models. The open triangle shows the 
position that T Tau would have if the water cooling were 
increased by a factor 3 and the OH cooling decreased by 
a factor 10. C-type shock models of Kaufman & Neufeld 
(1996) are shown in a grid as a function of density and 
shock velocity. J-type shock models (Hollenbach & Mc- 
Kee 1989) are shown as large dashed circles (see caption 
of Fig. 7). 



Because photodissociation of water by an UV field, 
that is not included in shock models, can convert water 
to OH, we suppose that the overabundance of the OH 
molecule is due only to the strong far-UV radiation field 
associated to T Tau (Herbig & Goodrich, 1986). The pho- 
todissociation cross Sect, of water at the Lya frequency is 
in fact ten times larger than the one of OH (van Dishoeck 
& Dalgarno 1984). A similar situation has been foimd in 
supernova remnants, where OH 1720 MHz emission is ex- 
plained as originated from C-type shocks, allowing that 
the action of an UV field creates sufficient OH from water 
dissociation (Wardle et al. 1998; Lockett et al. 1999). 

If this is the case, the far-infrared molecular emission 
spectrum of T Tau is primarily due to C-type shocks. To 
reconcile the OH observations with C-type shock models 
we need that the OH abundance, and thus its total cooling, 
were reduced by a factor 10 in favor of water cooling. From 
Table 6, if we reduce by a factor 10 the OH column density 
and we increase of the corresponding amount that one of 
the water, passing from 2-10"'^^ to 5.6-10^^, wc will increase 
the total water cooling by a factor of about 3. 
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Fig. 9. Total cooling from water lines versus total cooling 
from H2 lines, both normalized to the CO (high-J) cooling 
in T Tau and L1448 (Nisini et al. 1999b) (filled triangles) 
and according to shock models. The open triangle shows 
the position that T Tau would have if the water cooling 
were increased by a factor 3. C-type shock models of Kauf- 
man & Neufeld (1996) are shown in a grid where density 
increases upward and shock velocity from the left to the 
right. J-type shock models (Hollenbach & McKee 1989) 
are shown as large dashed circles (see caption of Fig. 7). 

Increasing the water cooling by this factor and de- 
creasing the OH cooling by a factor 10 would move the 
position of T Tau in the three diagrams of Fig. 7, 8 and 

9 in positions fully consistent with C-type shock mod- 
els (see the open triangles in the figures), at pre-shock 
densities of about l0^-^~^ cm~^ and shock velocities of 

10 < Vg < 20 km s"-'^. Comparing the pre-shock densities 
of these models with the densities derived from our LVG 
models, we found that a moderate compression factor (< 
10) would be required. In support of the fact that the 
molecular emission is due to C-type shocks is the evidence 
of the presence of strong magnetic fields in the outfiow re- 
gion associated with T Tau S (Ray et al. 1997). 

Our conclusion that the observed far-infrared molec- 
ular emission from T Tau can be explained by C-type 
shocks and that the atomic emission is probably origi- 
nated in J-type shocks is in agreement with the findings 
of van den Ancker et al. (1999). Their models, based on 
the near to mid-infrared II2 emission imply two temper- 
ature components at 440 K and 1500 K, which again are 
roughly in agreement with our finding of two components: 
one ranging from 300 K to 900 K and another at an higher 
temperature traced by the higher 3^p CO transitions. 



As to the origin of the C-typc shocks, responsible for 
the observed far-infrared molecular emission, we know 
from the LVG models that the emission region has a size 
of only few hundreds AU, assuming a spherical geometry. 
This implies that the shocks occur in a very compact re- 
gion, presumably very close to the binary system. There 
arc at least three mechanisms not mutually exclusive to 
explain the origin of the shocks: 

1. from the interaction region of winds coming from the 
two stars; 

2. from disk accretion on the youngest component of the 
binary system; 

3. from the interaction of the stellar wind with the molec- 
ular material in the circumstellar envelope. 

The first possibility seems the best one for originating fast 
(vg ~ 50 km s~^) dissociative J-type shocks, because the 
wind interaction would occur close to the stars where ve- 
locities arc siipposed to be high and the stellar field strong 
to dissociate molecules. The second possibility has already 
been suggested by Quirrenbach & Zinnecker (1997) to ex- 
plain the near-infrared H2 extended emission. The third 
one is probably at work in any case, because it is needed 
to explain the strong emission from CO, H2O, OH, as well 
as that from H2 (van den Ancker et al. 1999). 

5. Conclusions 

To summarise our results, we list the main findings of this 
study: 

1. The far-infrared spectrum associated to the binary 
system of T Tau shows strong emission from CO, 
H2O and OH molecules. 

2. Optically thin CO emission lines from high-J transi- 
tions are used to constrain the physical regimes of the 
gas: T-300-900 K and nn^ -IQ-'^^^cm--'^. The detec- 
tion of CO lines with Jup of 24 and 25 seems to indicate 
that a warmer component is also needed. 

3. H2O and OH emission is consistent with such condi- 
tions, however the higher excitation lines at the shorter 
wavelenghs are not well fitted by these models, indicat- 
ing that higher temperature gas should also be present, 
in agreement with the CO emission. 

4. Prom the assumption that all the far-infrared molecu- 
lar emission observed originate from the same region 
a very compact size of 300 - 400 AU of diameter is 
derived. 

5. The detection of [CII]158/im off-source at large dis- 
tances from T Tau shows that most of the ionised car- 
bon emission (^ 80%) is from an extended region and 
not originated in the vicinity of T Tau. The intrinsic 
emission in the LWS beam centered on T Tau implies a 
ratio [OI]63/xm/[CII]158/um of 115, much in excess of 
what is expected from photodissociation region mod- 
els. This emission is probably due to J-type shocks. 
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however it is not ruled out the possibility of a contam- 
ination by a PDR. 

6. The total cooling observed in the various species is 
comparable to the mechanical luminosity of the out- 
flow, indicating that the stellar winds could be the ulti- 
mate responsible of the line excitation through shocks. 

7. The study of the cooling from the main radiative com- 
ponents ([01], CO, H2O, OH and H2) is used to show 
that our data are consistent with C-type shock mod- 
els, only if we assume an overabundance of the OH 
molecule (by a factor 10) probably due to the intense 
far-ultraviolet radiation field. 

8. We finally consider the origin of the shocks responsi- 
ble for the far-infrared emission: colliding winds from 
the components of the binary system, disk accretion 
in the more embedded star and wind interaction with 
the molecular ambient cloud. Due to the large beam 
size of the LWS instrument, we are unable to specify 
which of the processes is the dominant one. 
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